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Abstract 

We present studies of a simple Anderson model Hamiltonian for Ce 3+ ions 



■ in cubic symmetry with three configurations (f°,f ,/ 2 )- Our Hamiltonian 
ON, 

includes: (i) One-channel Anderson model; (ii) Two-channel Anderson model. 

00 

■ Using the third order scaling analysis, we study stability of the non- Fermi 

10 ■ ■ , ■ , - 

ON . liquid fixed point of the two-channel Kondo model for Ce 3+ ions in cubic 



symmetry against the one-channel Kondo interaction. Using the non-crossing 



a 

i-Q ■ approximation (NCA), we also report detailed studies of our simplified model 

of one-channel & two-channel Anderson model which exhibits competition 



between the Fermi-liquid fixed point of the one-channel Kondo model and the 
non-Fermi fixed point of the two-channel Kondo model. We provide the phase 
diagram in the model parameter space and study the thermodynamics and the 
transport properties of our simplified model Hamiltonian. Thermodynamics 
and transport coefficients show the distinct behaviors for different numbers of 
channels. We confirm in detail that the NCA is a valid numerical method for 
the overcompensated multichannel Sj = 1/2 Anderson models. Our model 
study might be relevant to the non- Fermi liquid alloy Cei_ x La :r Cu2.2Si2. 
PACS Nos. 74.70.Vy, 74.65. +n, 74.70.Tx 
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I. INTRODUCTION 



The Kondo effect [|T| has been of great interest in condensed matter physics since its ob- 
servation. The proposed model Hamiltonian, a magnetic Si = 1/2 local moment interacting 
with the conduction electron gas, looked very simple but was non-trivial due to many body 
nature of the problem. Ever since, many generalized models have been studied to extend our 
understanding and to relate to real materials. The simplest Si = 1/2 orbitally nondegenerate 
Anderson model and s-d exchange model are now well understood for a single impurity 
case using several techniques. The numerical renormalization group (NRG) B method was 
able to provide complete information about crossover from the high temperature fixed point 
to the low temperature fermi liquid fixed point for these models. Subsequently, the exact 
diagonalization of these models was realized by the Bethe Ansatz (BA) Q, which also gives 
the exact solution for thermodynamics of these models. However, it has not proven possible 
to compute dynamical properties with the BA. Through non-crossing approximation (NCA) 
||, dynamics as well as thermodynamics |J have been extensively studied for the infinite 
on-site Coulomb interaction models. The Quantum monte Carlo method (QMC) has 
also been applied to study statics and dynamics for the simple Si = 1/2 models. Recently 
conformal field theory (CFT) |S||| has been used to study all properties near the low tem- 
perature fixed points. However we are still far from a complete understanding for realistic 
models which, for example, include the strong spin-orbit coupling, the crystal electric field 
effects, and multiple (more than two) configurations. In this paper, we study a realistic 
extension of the conventional simple approach to the Kondo effect for Ce 3+ ions, including 
the strong spin-orbit coupling, the crystal electric field effect, and multiple configurations in 
their simplest form. 



Ever since Nozieres and Blandin [|1(J introduced the multichannel Kondo model, its re- 
alization in real materials has been controversial. On the theoretical side, the multichannel 
Kondo models are well understood [^.|TT|-p^| irrespective of the experimental situations. 
When the degenerate channel number (iV c h), which is the number of the conduction electron 
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"bands" coupled to the impurity site, is greater than 2Sj (the impurity spin), the impu- 
rity magnetic moment is overscreened and the strong coupling fixed point (Fermi liquid) 
becomes unstable leading to a non-trivial fixed point (non- Fermi liquid behavior). In the 
underscreened and the completely screened cases, N c ^ < 2Si, the strong coupling fixed 
point is stable resulting in a Fermi liquid ground state. The spin susceptibility xCO an d 
specific heat coefficient C(T)/T for the two-channel Si = 1/2 magnetic Kondo model are 
proportional to log(2V/T) at low temperatures P,P,[iTl-|l3|1 , where Tk is the Kondo energy 
scale. The dynamic susceptibility shows Marginal Fermi liquid behavior ||14|| . The resistivity 
increases logarithmically and saturates to a constant with p(T) = p(0) [ 1 — cl^Jt/Tk ] with 
decreasing temperature ||. On the other hand, the one-channel Si = 1/2 Kondo model 
leads to the Fermi liquid ground state. In that case, magnetic susceptibility x(T) and spe- 
cific heat coefficient C{T)/T saturate to constants of order l/T K ||f|]. Resistivity increases 
logarithmically and saturates to a constant with p(T) = p(0) [ 1 — a(T/Tk) 2 ] with decreasing 
temperature []I5] . 

In this paper we study a model Hamiltonian for Ce 3+ ions in cubic metals with three 
configurations (/°, / , f 2 ). The nominal ground configuration f 1 can fluctuate to f° and f 2 
configurations by hybridizing with the conduction electrons. A one-channel Anderson hy- 
bridization interaction is present between f° and f 1 configurations. A two-channel Anderson 
hybridization interaction is present between f 1 and f 2 configurations. We report detailed 
studies of our simplified Hamiltonian of one-channel & two-channel Anderson model using 
the NCA. This simple model is quite intriguing in that we can study the competition be- 
tween the two different Kondo effects, that is, the Fermi liquid and the non-Fermi liquid 
fixed points. The distinct ground state physics for the different numbers of channels is clas- 
sified using the zero temperature analysis of NCA integral equations and the third order 
scaling theory. We calculate the thermodynamics and the dynamics of this simple model 
and find that all the calculated physical quantities show the behavior appropriate for the 
accessible different channel numbers. The static magnetic susceptibility displays a scaling 
behavior agreeing with the exact Bethe ansatz results in the two- and three-channel cases. 
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NCA calculation of entropy and specific heat is also compared with the Bethe ansatz results. 
The resistivity shows the correct temperature dependence near zero temperature agreeing 
with the conformal field theory results in the two- and three-channel cases. The sign and the 
magnitude of the thermopower are dependent sensitively on the relevant channel numbers. 
The peak position in the dynamic magnetic susceptibility is almost linear in temperature 
in the overscreened cases. A short paper which presents some of these results will appear 
elsewhere Jl6| . 

Our study is motivated by a recent discovery of non-Fermi liquid system, 
Ce x Lai_a;Cu2.2Si2 (x = 0.1) Hl7| . We summarize the experimental findings of this alloy 
system. The logarithmic divergence in both the magnetic susceptibility x{T) and the spe- 
cific heat linear coefficient 7(T) have been observed for Ce.rLa1_.jCu2.2Si2 (% = 0.1). The 



two-channel Si = 1/2 magnetic Kondo physics P,P,pT]-|T^] provides a theoretical framework 



to explain the thermodynamics of this system at low temperatures. j(T) initially increases 
in the presence of the magnetic field, which qualitatively agrees with the two-channel Kondo 
effect coming from the lifting of residual entropy [IBf. In the one-channel Kondo effect, the 



Sommerfeld coefficient decreases in the magnetic field due to the destruction of the Kondo 
effect. The Wilson ratio is estimated to be R m 2.7 from the slopes of two curves 
and j(T)), which compares well with theoretical value 8/3 for the two-channel magnetic 
Si = 1/2 Kondo model ||. The good agreement between the theoretical and the experi- 
mental Wilson ratios supports our crystalline electric field energy scheme described below. 
This system is pseudo-cubic (i.e., the crystal field scheme on the Ce 3+ site appears cubic). 
The best superconducting system with excess Cu shows almost isotropic magnetic suscep- 
tibility flSfl . That the pseudo-cubic T 7 magnetic doublet in J 1 lies lowest is inferred from 



the neutron scattering experiments ||19|| . The thermopower for CeCu2Si2 changes sign near 
70 K and stays negative below with a large extremum (—20 to — 30/xV/K) ||18|| , suggesting 
the presence of strong hole resonance scattering. As we will show below, these thermopower 
results also support our interpretation of the two channel magnetic Kondo physics. Though 
other experiments (e.g., specific heat and magnetic susceptibility) support the interpretation 



of them in terms of the two-channel Kondo effect, a linear temperature dependence in the 
resistivity remains as a puzzle. The \/T behavior in the resistivity is predicted from con- 
formal field theory treatment of the two-channel Kondo models ||,|9[]. In Fig. [I], we present 
our numerical calculation of resistivity and the experimental results measured in the alloy 
system Ce x Lai_ z Cu 2 .2Si2 (x = 0.1) |T7[ . It can be seen that the data curves downwards 
at lower temperatures which may indicate a crossover to a new fixed point. From all these 
experimental findings, we believe that the alloy system Ce x Lai_ x Cu2.2Si2 (x = 0.1) is a 
strong candidate for two-channel Si = 1/2 magnetic Kondo system. 

In addition, the two-channel, Si = 1/2 Kondo effect may be realized in other materials: 



two-channel quadrupolar Kondo effect [20| in some U alloys and two-level systems [21 



in 



metallic point contacts. U alloy systems include U0.2Y0.sPd3 p2| , U x Thi_ x .Ru2Si2 p3 



UCu 3 . 5 Pdi.5 1H, U .iTho.9Ni 2 Al3 Uo.iPro.gNisAla |§, U x S Cl _ x .Pd 3 @, Uo. 9 Tho.iBe 13 
27], and lX r Thi_. r Pd2Al3 | 29|| . All the above systems show a logarithmic divergence at 



low temperature in the linear specific heat coefficient and a square root saturation in the 
static magnetic susceptibility. U x Th 1 _ a .Ru2Si2 [2~3] is an exception showing a logarithmically 
divergent magnetic susceptibility. The two channel quadrupolar Kondo physics has been 
invoked to explain non-Fermi liquid behavior in thermodynamic and transport properties of 
U^Yi-^Pds for x = 0.2 p2j and other U alloy systems. Recently, the resistivity in a metallic 
constriction was observed to obey \/T behavior and was interpreted as due to two-channel 
Kondo scattering from atomic two-level tunneling systems |30"|j3T|| . 

Our paper is organized as follows. In section II, we introduce our simple model Hamil- 
tonian and analyze this model using the third order scaling equation. We briefly introduce 
the NCA in section III. Zero temperature analysis of NCA integral equations follows in 
section IV. In section V, we present the detailed numerical analysis of our simple model 
Hamiltonian. We conclude in the section VI. 
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II. MODEL HAMILTON! AN 



The single impurity Anderson model has been very successful || in describing Kondo 
systems (meaning magnetic transition metal elements embedded in normal metals, and 
dilute rare earth or actinide alloys). The thermodynamics is rather well explained by the 
single impurity properties for even highly concentrated Ce alloys PJ?2"|j. Coherence effects, 
arising from the lattice of Anderson or Kondo ions at low temperatures, do not play such 
an important role in thermodynamics. Transport properties also are well explained by 
the single impurity model except for the low temperature regime where a coherence effect 
prevails due to the coherent Bloch state formation leading to vanishing resistivity at zero 
temperature (the residual resistivity is larger than the room temperature one in the dilute 
impurity limit). This one-channel Anderson model can explain the complete screening of 
the magnetic moment at the local moment sites leading to the local Fermi liquid ground 
state discussed in the introduction. 

The two-channel orbital Kondo model (quadrupolar or two-level system) keeps the chan- 
nel symmetry guaranteed by time reversal symmetry, but can suffer from the ground state 
degeneracy lifting due to the Jahn- Teller effect. In general, the orbital Kondo model has the 
exchange anisotropy. It has been shown that the exchange anisotropy is irrelevant for the 
two-channel Sj = 1/2 Kondo models As will be shown in other publication [33|, the 
two-channel magnetic Kondo model for Ce 3+ is vulnerable to the channel asymmetry due 



to its orbital nature in the channel degrees of freedom ||10|| . However, the channel-mixing 
Kondo interaction can save the two-channel physics [[34] . 



The conventional theory of Ce 3+ impurities assumes the infinite on-site Coulomb inter- 
action, which removes the f 2 configuration from consideration 0, and as a result has no 
chance to get the two-channel Kondo effect which we shall explain below. When we relax 
our assumption about the strong on-site Coulomb interaction and we include the detailed 



atomic energy structure, we can develop a variety of model Hamiltonians [34 



In our simple model, we assume the magnetic f 1 J = 5/2 T 7 CEF doublet lies lowest in 
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the f 1 configuration, and we keep the two excited states - a singlet f° and the nonmagnetic 
f 2 T 3 CEF doublet. We find the one-channel Anderson model in mixing between f° and f 1 
configurations and the two-channel Anderson model in mixing between f l and f 2 configu- 



rations |35| . Other interesting Kondo interactions JJJ] arise when the excited triplets in the 
f 2 configuration are included. According to the group theory analysis, the hybridization is 
mediated only by the cubic T 8 conduction electrons between f 1 and f 2 (T 3 <g> T 7 = r 8 ), and 
by Ti between f° and f 1 for the mixing potential allowed in the cubic crystal. CEF states 
are schematically drawn in the Fig. ^| for this simple model. To see the essential physics, we 
restrict our attention to the simple case of isotropic hybridization and a free conduction band 
with Lorentzian/ Gaussian density of states (DOS). In this simple case, two components of 
conduction partial wave j c = 5/2, 7/2 can mix with the atomic orbitals. For the moment, 
we will consider only one angular momentum component, say, j c = 5/2 of the conduction 
band for our model study. We shall examine the effects of relaxing this in the other paper 



34| . Our model Hamiltonian is 

H — H ch + H at + Hi, (1) 

#cb = £ €C lr 7 a C eT T a + £ eC lr a na C eT 8 na, (2) 
ea enct 

H at = e \f ></°| + eiE l/ lr r« >< + ^ £ If^n >< f 2 T 3 n\, (3) 

a n 

Hi = Hi £4 7 J/° >< f l T 7 a\ + h.c. 

ea 

+Vn £(-ir +1/2 4r 8 J/^ra >< / 2 r 3 n| + h.c. (4) 

ena 

c eT 7 a{c\v 7 a) an d c er 8 na(4r g na) are the annihilation (creation) operators for conduction elec- 
trons of T 7 and r 8 , respectively, a =|, { denotes the time reversal pair for the irreducible 
r 7 representation, and n = ± is the quadrupolar index for T 3 and here acting as the chan- 
nel indices. The conduction electrons are assumed to be described by the uncorrelated 
Lorentzian/ Gaussian density of states with bandwidth D = 3 eV. eo,i,2 are the configuration 
energies for the empty (e = 0), single, and double occupancies, respectively. We lumped 
all Clebsch-Gordon coefficients into the hybridization constants, Voi and V12, except for the 
phase dependence on the cubic degeneracy indices. The phase dependence on the Kramers 



doublet index comes from time reversal symmetry. We will consider these two hybridization 
constants to be independent of each other to probe the competition between the one-channel 
and the two-channel Kondo physics. 

When the real charge fluctuations are removed from the model system in the Kondo 
limit, we have to construct the tensor operators representing each CEF states for the f 1 
configuration and the projected conduction electron CEF states. In this paper, the relevant 
tensor operators are for r 7j8 CEF states. We can show using the character table that 

r 7 ®r 7 = r 1 ©r 4 , (5) 
r 8 ® r 8 = Tt © r 2 © r 3 © 2r 4 © 2r 5 . (6) 

In the direct product, the first CEF states are written as ket, and the second as bra. The 
T 7 tensor operator (2x2 tensor) is a direct sum of charge operator (ri) and spin operator 
(r 4 ). Indeed, the Schrieffer- Wolff transformation leads to two interaction terms: the spin 
exchange interaction and the pure potential scattering term. The relevant terms of r 8 tensor 
operators are 2r 4 in our model. In the conduction electron T 8 tensor space, one of the two 
T 4 's gives rise to the ordinary S c = 1/2 spin operators with two degenerate orbital channels 
and the other to the S c = 3/2 spin operator with one channel . 



Hence, there are three distinct channel labels for conduction electron partial wave states 
about the Ce impurity in this simple model. One channel is just the T 7 doublet. The other 
two are the T 3 ± "orbital" states of the T 8 quartet. Each T 3 ± orbital has a T 7 "spin" doublet 
(recall r 8 = T 7 ©r 3 .) As shown in Fig. ^ for J = 5/2 conduction partial waves, the + orbital 
is "stretched" along the quantized axis (one of the three principal cubic axes, taken to be 
z here for definiteness. ) The "— " orbital is "squashed" in the xy plane. We note that the 
simplest example of T 8 partial wave quartet is, for zero spin-orbit coupling, ci-wave states 
with H > 3z 2 — r 2 and > x 2 — y 2 . The "spin" index is then real spin of the electrons. 

In the Kondo limit with the stable f 1 configuration, we may remove the f° — f , f — f 2 
charge fluctuations from the Hamiltonian of Eq. (4) using the Schrieffer- Wolff transformation 



1 37| to find the effective Hamiltonian. 



H x = JiS cV7 {Qi) ■ S F7 + ^r 8 n(0) • S F7 (7) 



n=± 



Ji = W 4 = W (8 ) 

-ei e 2 - ei 

^v = ^EI/ 1 ; r r«>^</ 1 ;r 7 ^|. (9) 

a/3 

5r 7 is the J 1 pseudo spin. S c r 7 (0) and S c r s ±(0) are the conduction electron pseudo spin 
densities at the impurity site of symmetry Tj and Tg, respectively. When the f 1 configu- 
ration is stable, its pseudo spin is coupled to the conduction band in a one-channel via f° 
configuration, and in a two-channel via f 2 configuration. Their couplings are both antiferro- 
magnetic. The unique feature of our Hamiltonian is that it can generate the one-, two-, and 
three-channel ground states depending on the model parameters. The competition between 
the Fermi liquid fixed point and the non-Fermi liquid fixed point can thus be investigated 
using this model Hamiltonian. 

We first analyze our simple model Hamiltonian using the third order scaling argument, 
i.e., the perturbative renormalization group (RG) [fLCjj3"8ll . At temperature T, only the 
conduction electrons (thermally excited) inside the band of order T with respect to the Fermi 
level play an important role in determining physical properties. Thus we can integrate out 
the band edge states (virtually excited states) to find the effective Hamiltonian. Though the 
following analysis is restricted to uj <C D(the conduction band width) and the perturbative 
regime (weak coupling limit), we can derive qualitative results out of this. For quantitative 
results, a full numerical renormalization group (NRG) study is required. 

It can be deduced from the scaling theory that the low temperature and the low energy 
physics is dominated by the one-channel or the two-channel Kondo effects depending on their 
relative magnitude of the antiferromagnetic couplings. We introduce an exchange coupling 
matrix in the channel space which is convenient for the derivation of the scaling equations. 
We can thus rewrite the one-channel and two-channel Kondo model in the following form 
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Hi = J ® S c (0) ■ S h (10) 
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/ J i \ 



(11) 



J 2 

V o o j 2 j 

Here S c and Sj are 5 = 1/2 operators. The scaling equations of our simple model Hamilto- 
nian up to the third order diagrams of Fig. [| are 

g = g 2 -^gTr[g 2 ], (12) 
g = iV(0)J. (13) 

The scaling equations in components are 

^ = g\ - \gi [ g\ + 2gi ], m = N{0) J x > o, (w) 
§| = gl - I 92 [ g\ + 2<?2 1 92 = N(0)J 2 > o. (15) 

Here x = log(D/T). We can identify three fixed points related to one-, two-, and three- 
channel Kondo physics. The one-channel, strong coupling fixed point {g*,g^) — (oo,0) is 



stable leading to the Fermi liquid ground state p8 |. The three-channel fixed point (2/3, 2/3) 
is stable along the line g± = g 2 in the g\ — g 2 plane, but unstable for any small perturbation 
from gi = g 2 . Finally, the two-channel fixed point (0, 1) is stable leading to the logarithmi- 
cally divergent thermodynamic properties at zero temperature. From the scaling analysis, 
we can infer the ground state physics: one-channel for J\ > J 2 ; two-channel for J x < J 2 ; and 
three-channel for J\ = J 2 . As will be shown in section IV, the zero temperature analysis of 
the NCA equations leads to the same conclusion. 

We now discuss the neglected other 8 T 3 irreps in the f 2 configuration [Q. The 9 T 3 
CEF states all contribute to the enhancement of the two-channel exchange coupling between 
/ 1 T 7 spin and the r 8 conduction electron spins. 

Hi = EE(-l) Q+1/2 ^ 2 4v sna \f^ra >< fr 3 n\ + h.c. (16) 

e ina 

The NCA can treat this problem with the extension that now the f 2 T 3 Green's function 
becomes an 9 x 9 matrix. See the Appendix for more details. The Schrieffer- Wolff transfor- 
mation leads to 
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(17) 




(18) 



Here e| is the energy level for the z-th / 2 r3 state. Hence multiple 1^ states in the f 2 
configuration lead to the enhancement of the two-channel exchange coupling. 

Particle-hole asymmetry in the conduction band density of states (DOS) is also important 
in determining the ground state weights of the one-channel (/° — f 1 ) and the two-channel 
(f 1 — f 2 ) contributions. In the scaling approach, the particle-hole asymmetry is completely 
neglected. However the NCA can take into account this particle-hole asymmetry. The 
occupied conduction electron states (hole) contribute to the f° self energy (see the section 
III), while the unoccupied states (particle) to the f 2 self energy. Hence the more weight 
in the particle side can enhance the effective hybridization strength between f 1 and f 2 
configurations. The particle-dominant conduction band DOS will lead to the enhancement 
of the two-channel exchange coupling. 



system - one-channel & two-channel Anderson model. Though this model is highly simplified 
compared with the full model p4| , we can simulate the full model using this simple model 
Hamiltonian to study the different physics for varying channel number. 

In the NCA, our starting basis is the conduction band plus the atomic Hamiltonian pro- 
jected to the atomic electron Fock space and treat the hybridization between the conduction 
band and the atomic orbital as a perturbation. The strength of this approach is that the 
strong on-site Coulomb interaction for atomic electrons is treated accurately at the outset. 
Pseudo particle Green's functions are introduced for each atomic electron occupation state 
which is neither fermionic nor bosonic. The price we pay is that we cannot apply the con- 
ventional Feynmann diagram techniques to this strongly correlated problems. Thus special 



III. NON-CROSSING APPROXIMATION (NCA). 



We now apply the non-crossing approximation (NCA) |J to study our simple model 
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Green's function technique was developed by many investigators P]^,j39H42"f ■ This approach 
may be justified by expansion in 1/N which reorders the diagrams by treating NV 2 as 0(1). 
Here N is the nominal atomic ground state degeneracy and V is the hybridization strength 
between the conduction band and the atomic orbitals. In the NCA, pseudo particle self en- 
ergy diagrams include only the leading order skeleton (non-crossing) diagrams and they are 
solved self-consistently. For the one-channel models, this theory includes all the diagrams 
up to 0(1/N) order and a subset of higher order diagrams up to the infinite order. Though 
vertex corrections, which is 0(1/N 2 ), are not included, this approach is a conserving ap- 
proximation |40|| . For the overscreened multi-channel Anderson models, it has been shown 
|14|| that the 1/N approach becomes exact in the limit M,N — > oo with fixed M/N ratio 
(M is the number of channels). 

When we study the most general three-configuration model, the same symmetry con- 
duction electron can be involved in the two mixing processes, e.g., f° — f 1 and f 1 — f 2 for 
Ce 3+ atoms. Generally, a specific vertex correction is required to get the right Kondo energy 
scales in this case. Recently, vertex corrections were included in the study of the finite U 
spin 1/2 Anderson model P3| . A simplifying feature of our model Hamiltonian is that the 
leading vertex correction vanishes, since two different symmetry conduction electrons are 
involved in the hybridizations f° — f 1 (T?) and f 1 — f 2 (Tg). This feature greatly simplifies 
the numerical work and formalism. 

From the leading order skeleton diagrams of Fig. [|, we find the self-consistent integral 
equations: 



EoC*) = — E /<*e N{e)f{e)G 1 {z + e) : (19) 

TV a J 

Ei(z) = ^ fdeN(e)f(-e)G (z-e) 

71 J 

£ / deN(e)f(e)G 2 (z + e), (20) 

= — E / rfe #(e)/H0Gi(* - e), (21) 

G N (z) = VtT ? r 4 ^vriV(0)|^| 2 . (22) 

Z — Cat — L,n{ Z) 



£12 

7T 
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Here N(e) is the normalized conduction band DOS at the fermi level such that N(0) = 1. 
So,i,2(- 2 ) and Go,i,2(z) are the self energy equations and Green's functions for /°, f 1 ^, and 
/ 2 T 3 atomic states, respectively, /(e) is the Fermi-Dirac distribution function. is the 
hybridization strength characterizing the width of the renormalized atomic electron spectral 
function peak. One of the strong points of the NCA approach is that we can easily study 
any form of the conduction band DOS, as opposed to other approaches, e.g., Bethe ansatz 
or conformal field theory, to Anderson/Kondo models. As a concrete example, we will use 
the structureless Lorentzian/ Gaussian DOS to see the many body physics. We solved these 
coupled integral equations numerically to study the thermodynamics and the dynamics of 
the model Hamiltonian. 

Pseudo particle Green's functions are not directly measurable. All the physically mea- 
surable quantities are given by the convolution of the pseudo particle Green's functions. 
Now it is convenient to introduce the spectral function (An(u)) for each pseudo particle 
Green's function and its negative frequency spectral function (a N (u)). 

A n (uj) = --Im 1 (23) 

a N (co) ee e-P" A N (co). (24) 

Since ajv(c<j) always appears in combination with the impurity partition function Zf, there 
is arbitraryness in overall factor in its definition. The above NCA integral equations do 
not have a N {uj) in them. Whenever we calculate any measurable physical quantities, the 
spectral functions appear with the Boltzman thermal factor divided by the impurity partition 
function. In the numerical work, we calculate aAr(u;)'s self-consistently. Roughly speaking, 
ajv(c<j)'s are occupancy distribution function weighted by the thermal factor. 
The impurity partition function for our simple model is defined by 

-Zf ee J du [a (u>) + 2ai(u) + 2a 2 (w)] . (25) 

This partition function includes the many body effect of the interaction between impurity 
and the conduction band and is exact in its form. The approximation comes in when we 
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choose the approximate Green's functions for the atomic states. 

In our simple model, only two symmetry electrons are present: T 7 and T 8 . Their spectral 
functions (measurable) are defined as a convolution of two neighboring configurations' spec- 
tral functions. It can be shown that these two interconfiguration atomic spectral functions 
are given by 

Pr» = - / d( a (C) Ax(C + ") = PoiM, (26) 

p r » = 1 + e z - J d( ai(C) A 2 (( + uj) = paiu) (27) 



in this conserving approximation |)4T) |. From now on, we will use the notations poi, pu in 
favor of pr r ,Pv 8 - This approximation does not include any vertex correction. 

From the leading bubble diagram || , the static magnetic susceptibility per Ce 3+ impurity 

is 

1 75 
X (T) = - p 2 cS Xf(T), & = — (28) 

X{T) = ~Jdt ai(C,r) ReG^CT). (29) 

Here /ib is Bohr magneton. The reduced dynamic magnetic susceptibility is 

» 1 — e _/3a; r 

X (w, T) = - y rfC ai(C, T) A,(C + u,T). (30) 

We are going to calculate these physical quantities and transport coefficients. 

IV. ZERO TEMPERATURE ANALYSIS. 



Zero temperature analysis fl4]j44}| of the NCA integral equations leads to a qualitative 



understanding of our model system. Recently, conformal field theory approach || calculated 
the scaling dimensions exactly for the multi-channel Kondo exchange models with any size 
of impurity spin. In Ref. [|I4]] , the evaluation of the scaling dimensions is extended to the 
multichannel Anderson/Coqblin-Schrieffer models using the functional integral formulation. 
Both models become congruent when the impurity spin is Sj = 1/2. 
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The main results of the zero temperature analysis of the NCA integral equations are: 
(1) We can find the criterion which predicts whether the ground state is that of the 1, 2, 
or 3 channel model. (2) The Kondo energy scale (To) can be estimated analytically for 
the case Toi = T^. To in the one- and two-channel model parameter regimes is shown to 
vanish as the f 2 T 3 energy level approaches that of the f° configuration. (3) We obtain the 
correct scaling dimensions for the overcompensated cases which agree with the conformal 
field theory results |§ . (4) The crossover physics between the parameter regimes for different 
numbers of channels can be identified. 

The self consistent NCA integral equations can be transformed into the differential equa- 
tions for the flat conduction band in the wide band limit: D 3> | ex,2 1 - We analyze the zero 
temperature NCA equations in the asymptotic limit \u — T | <C T . T is the threshold 
energy below which the pseudo particle Green's functions become purely real. We introduce 
the inverse Green's functions and transform the self energy equations at zero temperature 



into the coupled non-linear differential equations. [04 
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Then g and g 2 can be shown to be related by 



g 2 (-D) = D + e 2 . 
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In the zero temperature analysis, it is more convenient to define the "negative frequency" 
spectral functions by ai(u) = e _/3<JJ Ai(u)/Zf. These spectral functions vanish above the 
threshold energy Tq and satisfy 



4-[ a (^)|^oM| 2 ] = - — 

OLd 7T 
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ai(u), 
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ai(uj)\gi(uj)\ 2 ] = -— oq(u) - — ^ a 2 (w), (38) 
-f [a 2 M|£ 2 H| 2 ] = — ai(a;). (39) 

It can be shown from the above relations that 

[ a (uj)g (uj) + 2a 1 (u)g 1 (u) + 2a 2 (u)g 2 {uj) ] = a (u) + 2ai(a;) + 2a 2 (cu). (40) 

By integrating this equation from c<j = — oo to uj = E , we find the additional relation, 

[ a (u;)g (u;) + 2a 1 (u;)g 1 (u>) + 202(^)52(0;)]^^ = I. (41) 

This identity will be useful in finding the asymptotic behavior of the "negative frequency" 
spectral functions. As a corollary we have another identity in the asymptotic limit, 

aoMlfloM | 2 = a 2 Mlfl2M| 2 ^ 
Toi Ti2 

This relation can be proved by using the boundary condition at uj = E . We calculate 
the physical atomic spectral functions and the dynamic magnetic susceptibility defined in 
section IV. 

poiM = J de [ao(e)Ai(e + u) + A {e)a 1 (e + u)] , (43) 

P12M = / de [ai(e)A 2 (e + u) + A 1 (e)a 2 (e + u)] , (44) 

x"(w) = J de [ai(e)Ai(e + w) - Ai(e)ai(e + w)]. (45) 

From the dependence of the spectral functions near u = 0, we can infer the finite temper- 
ature dependence of transport coefficients as will be discussed below. 

We now discuss the phase diagram in the model parameter space using the zero tem- 
perature analysis. 7 C decides the low energy and the low temperature behaviors of our 
model Hamiltonian. 7 C measures the relative magnitude of the antiferromagnetic coupling 
strengths when the charge fluctuation is removed in the model Hamiltonian. Noting that 
E « ei + 0(V&, V? 2 ), we find 



7c 



ttN(0) 



1 1 

T 2 ~ Ji 



(46) 
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which illustrates the correspondence to the scaling analysis. If 7 C is greater than zero, the 
singular behavior shows up in the f° Green's function, and not in the f 2 Green's function. 
Hence the system will be dominated by the f° and f 1 sector leading to the one-channel 
Kondo effect. When 7 C is less than zero, the singular behavior shows up in the f 2 Green's 
function, and not in the /° Green's function. In this case f 1 and f 2 sector (the two-channel 
Kondo physics) determines the low temperature behavior of the system. When 7 C = 0, f° 
and f 2 become equivalent asymptotically — E \ <C T ). Both Green's functions develop 
the singular behaviors at the ground state energy. In this model parameter space, the 
three-channel Kondo model is realized. 

The characteristic Kondo energy scale T is found from an integration constant which 
connects the low and high energy states. We can obtain the integration constant for the 
case Ti2 = r i(= T). We will analyze this case in detail and indicate subsequently how to 
extend the zero temperature analysis to the general hybridization strength. 



A. The symmetric hybridization limit: Ti2 = roi(= T). 

When we take a symmetric hybridization limit T 12 = r i(= T), the equations are sim- 
plified. 

£2 



7c | , ■ 
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(47) 
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The last two relations hold true in the asymptotic limit, \oo — E \ <C T . The different ground 
states then are determined solely by the sign of the f 2 configuration energy relative to the 
f° configuration energy. Removing the variable cu, we can find the differential equations 
between the inverse Green's functions. Integrating from their values at uo = —D to those at 
u, we find the integration constant 
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in the wide conduction band limit, D » |ei i2 |. We identify three cases for evaluating T . 
(1) One-channel case: 7 C > or e 2 > 0. f° Green's function develops a singular behavior at 
the threshold energy, while f 2 Green's function does not. Thus the ground state physics is 
dominated by the sector f° — f 1 leading to the one-channel Kondo effect. The integration 
constant and the Kondo temperature are 
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Here A = T/n. We can find the asymptotic behavior of the Green's functions for each 
atomic state. 
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(54) 
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Though the scaling dimension is not correct, the estimated Kondo temperature is correct 
within order unity and x(0) ~ l/T . An interesting observation is that Kondo temperature 
vanishes as e 2 tends to zero, i.e., the three-channel parameter regime (see the discussion 
below). 

The detailed derivation of the asymptotic behavior is not included since NCA does not 
produce the Fermi liquid ground state in the one-channel model. We just give a brief 
summary of the zero temperature analysis which is relevant to our study. A 2 (uj) vanishes as 
Q 2/3 at the threshold energy, while A (Ai) is singular as fT 2/3 (iV 1/3 ) at uj = E . Thus, 
the physical spectral function p± 2 vanishes as \u\^ 3 at the Fermi energy, while poi is finite. 

(2) Two-channel case: 7 C < or e 2 < 0. In contrast to the one-channel case, the 
f 2 spectral function has a singular structure leading to two-channel ground state. The 
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integration constant and the Kondo temperature are 
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Note that the Kondo temperature vanishes with e 2 — > 0. We can find the asymptotic 
behavior following the standard zero temperature analysis. 



go{u) = g 2 {u) + |e 2 |, 



T n 



ft 1 / 2 [ 1 + dfi 1 / 2 + 0(6) ], 



^M = ftV2 [ 1 + c 2 n 1 ' 2 + 0(Cl) ], 



Eq — UJ 



Cl 



6 



2 1 + 



A 

M 



A 



c 2 = -- 



(60) 
(61) 

(62) 
(63) 

(64) 



The asymptotic behavior above the ground state energy E can be obtained from the ex- 
pressions below Eq by the analytic continuation. Furthermore, we find for the "negative 
frequency" spectra that 
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xi = -(2ci + c 2 ) 

X2 = g(Cl + 2C 2 ) 
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Note that £0,1,2 > 0. Now we can find the asymptotic behavior of the pseudo particle spectral 
functions. 
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As expected, A (u>) vanishes at the threshold energy and does not develop any singular 
behavior. On the other hand, Ai(u) and A 2 (uj) are singular at the threshold energy. Finally 
the physical spectral functions are given by in the asymptotic limit. 
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(77) 
(78) 
(79) 
(80) 
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Poi(u) vanishes at u — and increases as |a;| near the fermi level. This spectral depletion 
at the fermi level is also confirmed in the finite temperature NCA calculation. p 12 is peaked 
right at the fermi level and has more spectral weight below than above the fermi level since 
x 2 < x\ for To < A. The finite temperature NCA results also confirm this observation. 
The spectral functions become non-analytic at the Fermi level at zero temperature. Note 
that the dynamic magnetic susceptibility is step function like at u — 0, which is none other 
than marginal fermi liquid behavior [14| . From the asymptotic form of the zero temperature 
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spectral functions, we can infer the finite temperature dependence of the resistivity using 
the Kubo formula 



1 Tie 2 f ( df \ 



P(T) 

For u,T <C T , the relaxation rate for the conduction electron is approximately 
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Near the zero temperature, we may replace pi 2 (u,T) by our zero temperature one and 
find the y/T temperature dependence. In fact one complication arises due to the angular 
averaging. Since the angular harmonics conjugate to p i2 is positive definite, still the above 
simple argument applies. In the one-channel case, we cannot use this kind of simple argument 
to find the low temperature behavior of resistivity. The scaling dimensions agree with those 
from the conformal field theory approach. 

We summarize the low temperature properties: (i) The atomic spectral function p i2 is 
peaked right at the Fermi level independent of the occupancy, while p i vanishes at uj = and 
is depleted near the fermi level; (ii) The dynamic magnetic susceptibility is step function-like 



at zero frequency; (iii) The resistivity obeys the scaling behavior p(T) = p(0) [ 1 — a^T/T } 
near zero temperature. All these results will be shown consistent with numerical NCA 
calculations at finite temperatures. 

(3) Three-channel case: j c = or e 2 = 0. In this case, the f° and f 2 configurations are 
equivalent asymptotically leading to three-channel ground state. The integration constant 
and the Kondo temperature are 
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See the Appendix for the detailed derivation of asymptotic behavior. This model case is 
not different from the analysis of the general overcompensated models. We simply list the 
asymptotic behavior here. 
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The nagative spectral functions are 



1 



1 - x x - Vl ft 3 / 5 



ai (uj) = — fT 3 / 5 
01 

a 2 (cu) = ^ ft" 2 / 5 [l - x 2 tt 2 ' b - y 2 fl 3 / 5 + ■■■} 0(6), 

4 T 
xi = 2ci = -, yi= 0; x 2 = 0, 1/2 = 2rf 2 = 7^. 



And the pseudo particle spectral functions are 
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From the pseudo particle spectral functions, we can find the physical spectral functions in 
the asymptotic limit. 
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Here B(p, q) is the Beta function. The spectral functions become non-analytic at the Fermi 
level at zero temperature. One important observation is that the power laws agree with those 
obtained from the conformal field theory treatments for the overcompensated multichannel 
S — 1/2 models. Since poi and p i2 are equivalent in the asymptotic limit, the angular 
functions are factored out in the conduction electron scattering time. Thus we can read off 
the power law of the resistivity a = 2/5. 

Again, the low temperature properties are summarized below, (i) The atomic spectral 
function is peaked right at the Fermi level independent of the occupancy; (ii) The dynamic 
magnetic susceptibility diverges as uj~ 1 ^ at zero frequency; (hi) The resistivity obeys the 
scaling behavior p(T) = p(0) [ 1 — a(T/T ) 0A ] near zero temperature. All these results will 
be shown consistent with numerical NCA calculations at finite T. 



B. Crossover physics. 

We expect the smooth crossover from the high temperature regime to the low temperature 
regime where the relevant channel number physics shows up. We can estimate the crossover 
temperature in the one- and two-channel model parameter regime using the zero temperature 
analysis. Noting that the asymptotic analysis assumes go « e 2 for the one-channel regime 
and g 2 « 1 62 1 for the two-channel regime (see the integration constant), we can estimate 
the crossover energy scale below which one- or two-channel physics is realized. In the one- 
channel case, the integration constant can be rewritten as 
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We can see that the relative magnitude of go and €2 determines the differing channel behavior. 
The crossover energy scale can be defined by the relation g = e 2 . The crossover temperature 
is found in the one-channel case 
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And in the two-channel case, we find the crossover temperature 
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Substituting the Kondo energy scale, we can see that the crossover energy scales are 
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Note that the crossover temperature vanishes as | e 2 1 5 / 2 with e 2 — > 0. We intentionally write 
the crossover temperature as the prefactor times the corresponding three-channel Kondo 
temperature. If this prefactor is greater than 1, the system will not display the three-channel 
behavior with decreasing temperatures but will flow directly to the one- or two-channel fixed 
point. When this prefactor is less than 1, the system will display the three-channel behavior 
before finally flowing to each channel fixed point. In the one-channel regime, the system 
flows directly from the high temperature regime (local moment regime) to the one-channel 
fixed point when e 2 >> A with decreasing temperatures. In the opposite limit, e 2 << A, the 
system will flow from the local moment regime to the three-channel regime, and finally to the 
one-channel fixed point after passing through the crossover temperature. In the two-channel 
regime the same argument can be applied as in the one-channel case. 



C. General asymmetric hybridization case. 



Though the Kondo energy scale can not be estimated analytically for the general case 
Toi 7^ r i2 , still the asymptotic behavior (uj < E ) can be sorted out. Since the physics 
is the same as the above analysis, we shall just write the asymptotic properties of Green's 
functions. 
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(2) 7c < 0: 
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(3) 7 C = 0: In this case, f° and f 2 are equivalent asymptotically. 
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Though the Kondo energy scale cannot be evaluated analytically, it can be estimated from 
the numerical NCA calculation of, e.g., magnetic susceptibility and fitting it onto the uni- 
versal results. 



25 



V. NUMERICAL ANALYSIS. 



In this section we present results from our full numerical study at finite temperatures. 
We studied the model for the parameters listed in Table U2. This covers the one-, two-, and 
three-channel Kondo regimes. For simplicity, we have chosen the same hybridization strength 
in the f° — f 1 and the f l — f 2 sectors. The relevant channel-number-dependent physics can 
then be studied by varying the relative position of the f° and the f 2 configuration energies. 
This simple choice of the hybridization further makes it possible to find the characteristic 
Kondo energy scales analytically, which is described in the previous section in detail. 

Our main results are: (1) The magnetic susceptibility shows a scaling behavior and 
agrees well with the exact Bethe ansatz results in the two- and three-channel parameter 
regimes; (2) The NCA results of residual entropy in the two-, and three-channel models 
agree with the exact ones within order 0(1/N). The Kondo anomaly peak in the specific 
heat also agrees with the exact one in its magnitude; (3) The thermopower is a diagnostic 
to display the ground states for different numbers of relevant channels for our model; (4) 
The dynamic magnetic susceptibility is distinct between the one-channel model and over- 
compensated (two- and three-channel) models; (5) Due to the simplifying features of our 
model, the resistivity shows a bendover at low temperatures in the one-channel parameter 
regime. The resistivity in the two- and three-channel parameter regimes shows temperature 
dependences near T = OK agreeing with the conformal field theory results; (6) We confirmed 
in detail that NCA is a valid numerical self-consistent non-perturbative method in studying 
the overcompensated multi-channel Sj = 1/2 Anderson model at T > 0. 



A. Entropy and specific heat. 

The entropy and the specific heat due to the magnetic impurity can be calculated from 
the free energy obtained in NCA through numerical differentiation. These thermodynamic 
quantities include very important information about the nature of the ground state. We can 
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estimate the characteristic energy scales in the Kondo models from the temperature variation 
of the entropy. In general, the entropy will increase with increasing temperatures, until the 
frozen impurity degrees of freedom are released. Our model Hamiltonian is expected to have 
the entropy S = /c#lri5 at high enough temperatures, while the residual entropy at zero 
temperature will show behavior expected for different numbers of relevant channels. 

The entropy and the specific heat are displayed in Fig. || and [7] for the model parameter 
sets studied here. In the one-channel case, the Kondo anomaly peak is well separated from 
the Schottky anomaly peak coming from the interconfiguration excitations. The Kondo 



anomaly peak has a magnitude comparing well with the exact results |13| for the Kondo 
exchange model. No residual entropy remains with T — > 0. In the two-channel case, the 
Kondo anomaly peak is not clearly separated from the charge fluctuation peak for most 



of our model parameters. The residual entropy agrees within 5 % with the exact one [13 
and the discrepancy can be explained by the 0(1/ N 2 ) and higher order corrections. To see 
this, we note that the entropy for spin 1/2 conduction electrons has an explicit dependence 
on the impurity spin degeneracy [§,0. Strictly speaking, the NCA results are valid for N- 
fold impurity spin and M-fold channel degeneracies as iV — > oo with N/M fixed. Hence 
it is natural to expect 0(1/ N 2 ) corrections to the entropy through the neglected vertex 
corrections. In the three-channel case, the Kondo anomaly peak is reduced further due to 
the increased residual entropy and is almost merged with the charge fluctuation contribution 
for model set 4. For model set 5, we can see a weak indication of separation. The residual 
entropy is increased compared to that of two-channel case and its magnitude is again a little 



bit smaller than the exact one [pl| due to the neglect of the higher order contribution in 
1/N expansion. 

B. Static magnetic susceptibility. 

The static magnetic susceptibility is a direct indicator of the nature of the ground state 
for the magnetic Kondo model. As is well documented, the magnetic susceptibility diverges 
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logarithmically x{T) oc log(To/T) for the two-channel Si = 1/2 magnetic Kondo model at 
zero temperature, and diverges algebraically for the three or the higher channel Si = 1/2 
Kondo model (x(T) oc [T /T] 1_2An ; A n = 2/(n+2): n is the channel number (> 3)) |ggg3 



The one-channel Kondo model has a constant magnetic susceptibility at zero temperature 
(x(0) ~ l/^o)- From the NCA, using the leading order bubble diagram 0, the magnetic 
susceptibility for one Ce 3+ impurity is given by 

1 75 

X (T) = -fa x(T), fa=^fa (121) 



X(T) = -|-/dC «i(C,T) ReG r 1 (C,T). (122) 



Here /ib is the Bohr magneton. 

Our numerical results for x(T) clearly show the right tendency for each possible low 
T channel number. The magnetic susceptibility in the one-channel regime, for the model 
parameter sets 6, 7, 8 shows approximate scaling behavior and clearly has the negative 
curvature at low temperatures indicating approach to the Fermi liquid ground state. The 
deviation of the scaling behavior at low temperatures seems to come from the pathological 
behavior of NCA in this one-channel case [Q . 

The x(T) curves in the two-channel regime, for the model parameter sets 1, 2, 3, also 
show scaling behavior and diverge logarithmically at low temperatures (Fig. ||). Our results 
are compared with those of Ref. ||13|| . The fitting to the exact Bethe ansatz numerical results 
13fl is quite good. We believe that the high temperature deviation comes from the M = 3 to 



M = 2 crossover physics described in the previous section. Note that Bethe ansatz results 
are for the pure two-channel s-d exchange model. To get the fitting to the Bethe ansatz, we 
slide the temperature axis to find T#- = 0.3 x T (TV from the Bethe Ansatz). Here T is 
the Kondo energy scale estimated from the zero temperature analysis. 

At high temperatures we cannot distinguish the different numbers of relevant channels 
clearly. Note that distinct physics for different numbers of relevant channels shows up at 
low temperatures below the crossover temperature which was estimated in the previous 
section. This observation is supported by results from the three-channel model parameter 
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sets. The high temperature deviation from scaling is very weak in this case which will not 
show the crossover physics. The magnetic susceptibility in the three-channel model also 
shows a scaling behavior. Since the three-channel case lies exactly on the boundary between 
the one-channel and the two-channel regime, we probed the three-channel case by varying 
the position of the bare f 1 configuration energy while the two excited configurations f°, f 2 
energies are kept equal. Fitting to the exact results is quite good. The Kondo energy scale 
To estimated from the zero temperature analysis agrees with that in the exact Bethe ansatz 
(TtfofRef EH). 



C. Spectral functions. 

The interconfiguration spectral functions show a distinct behavior for different numbers 
of relevant channels near the the Fermi level (u = 0) depending on the model parameters. In 
our simple model, we have two symmetry atomic spectral functions of T 7 and r 8 : pr 7 {&) = 
poi{u) and pr 8 M = Pvz{u). 
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p l2 (u) = J d( e-K A l (C) A 2 (C + u) (124) 



In the one-channel model parameter regime (see Fig. |9|), poi(^) develops the Kondo 
resonance peak just above the Fermi level and pi%{ui) is depleted near u = and tends to 
zero with T — > at u = 0. This confirms our zero temperature analysis in section IV. For 
comparison, note that the Kondo resonance amplitude is big compared to the two-channel 
or three-channel cases. 



In the two-channel regime (see Fig. |T^), pi2(^) is peaked below lu = and its peak 
position tends to uj = with decreasing temperature. Note that the Kondo resonance 
amplitude is reduced compared to the one-channel Kondo resonance amplitude. On the 
other hand, poi(^) is depleted near u = 0. 
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In the three-channel parameter regime (see Fig. [11]), the two spectral functions are equiv- 
alent asymptotically when uj — > — uj transformation is accounted for. Peak positions tend 
to zero with decreasing temperature. Both develop the Kondo resonance peak with further 
reduced amplitude. 

As mentioned above, the positions of the Kondo resonance peak (see Fig. |TJ) show a 
distinct behavior for different numbers of relevant channels. The peak position saturates to 
a constant value in the one-channel cases while it vanishes in the two- and three-channel 
cases with decreasing temperature. In addition, the detailed functional forms of the Kondo 
resonance peak are different depending on the relevant channel numbers. In the overscreened 
cases (M = 2,3), the peak structure becomes non-analytic with decreasing temperature as 
shown in the zero temperature analysis of the NCA integral equations. Finite temperature 
washes out this non-analytic behavior at the Kondo resonance peak. In the one-channel 
case, the atomic spectral functions remain analytic down to zero temperature. 

Since we are not considering all the atomic energy levels, the full atomic spectral functions 
(measured from the photoemission experiments) cannot be defined in our simple model. 
Though the high energy physics of real systems can not be properly treated within our simple 
model Hamiltonian, the low energy or low temperature properties can be studied using the 
restricted spectral functions. Note that low temperature/energy physics is governed by 
the Kondo resonance peak development. Though the spectral depletion is found in the 
interconfiguration spectral functions, we do not expect that the photoemission spectroscopy 
can observe this feature unless it can distinguish the atomic electron symmetry. Measurable 
atomic spectral functions are shown in Fig. [13| for each relevant channel case. 

D. Dynamic magnetic susceptibility. 

The dynamic magnetic susceptibility measures the magnetic excitation structure. Since 
the properties of the magnetic excitations are related to the interaction of the local magnetic 
moment with the conduction electrons, the relevant channel number will determine the 
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nature of the magnetic excitations. The dynamic magnetic susceptibility is expected to be 
distinct for each number of relevant channels. We have already seen this channel dependence 
in the static magnetic susceptibility. 

Dynamic magnetic susceptibility is defined as the spin-spin correlation function and can 
be measured directly from the neutron scattering experiments. From the leading bubble 
diagram, the reduced dynamic magnetic susceptibility is 



Our reduced static magnetic susceptibility is related to the above magnetic response function 
by the Kramer-Kronig relation: 



The neutron scattering experiments measure the structure function S(u, T) oc [ b{uj) + 
1 ] x"(c<j,T). The dynamic susceptibility can be quantitatively characterized by its 
'linewidth' dependence on temperature in addition to its overall functional shape. We may 
define the linewidth T(T) by the peak position of x'(u;,T), which can be measured directly 
in the inelastic neutron scattering experiments. 

The variation of x {u,T) and T(T) with temperature are displayed in Fig. |TJ| and [T5|. 
For the dynamic magnetic susceptibility, we display the reduced form x {u, T)/x {T[T), T) 
as a function of u/T(T). The distinct behavior for different number of relevant channels 
is clearly evidenced: in the one-channel regime, the x"(cu,T) curves converge with decreas- 
ing temperatures (not shown). This is clearly supported by the saturation of T(T) at low 
temperatures. In the two- and three-channel regimes, T(T) vanishes algebraically (close to 
linear) and the dynamic magnetic susceptibility does not converge (not shown) in contrast to 
the one-channel case, but instead develops non-analytic behavior at uj = 0. The reduced dy- 
namic magnetic susceptibilities (defined above, see Fig. [TJ]) roughly show a scaling behavior 
between two extrema. 

The physics of the linewidth is quite important in understanding the nature of the 
magnetic spin screening. The distinct behavior for different numbers of channels shows up at 




(125) 




(126) 
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low temperatures below To. The impurity spin flipping time (Yf) due to the hybridization will 
be given by the inverse of the linewidth, Tf ~ 1/T(T). On the other hand, thermally excited 
conduction electrons close to the impurity site will pass through it in a time r c of order of 
1/T from the uncertainty principle. At high temperatures above the Kondo temperature, 
Tf ^> t c . Thus the impurity spin rarely flips while the conduction electrons pass by the 
impurity site. Hence the Curie behavior is expected in the magnetic susceptibility. At low 
temperatures below the Kondo temperature, differing channel number physics shows up. 
In the one-channel regime, Tf <C r c and the impurity spin flips vigorously averaging out 
its spin to zero and leading to the Pauli behavior. In the two- and three-channel regimes, 
Tf pa t c and spin screening is not complete leading to the non-Fermi liquid ground state. 
This interpretation agrees with the low temperature behaviors of magnetic susceptibility. 



E. Transport coefficients. 

Using the Kubo formula in a dilute impurity limit, where the inter- impurity correlation 
can be neglected, we have calculated the transport coefficients, resistivity and thermopower. 
The anisotropic conduction electron scattering rate is 

r-'ikuj) = [ r 01 e^ 2 \k) Pol (u) + r 12 e 8 5/2) (^) Pl2 (u) ], (127) 

ef /2 \k) = ±- [ 6 - m ]; &f /2 \k) = -^-[Q + ], (128) 
$(fc) = 15 cos 4 9 - 10 cos 2 9 + 1 + 5 sin 4 9 cos 2 2<p. (129) 

Here n irav is the impurity concentration. Note that the conduction electron scattering rate 
contains the hybridization since the interconfiguration spectral functions are involved. The 
crystal harmonics are normalized such that the integration of ®y (k) and Bg 5 ''' 2 ^) over the 
solid angle leads to 1 and 2, respectively. Note that the the crystal harmonic Q^ 2 \k) van- 
ishes at "hot spot" angles of (0, (p) = (0,-), (tt,-), (tt/2,0), (vr/2, vr/2),(7r/2, ic),(ir/2, 3vr/2), 
while Q^ 2 \k) is positive definite. This feature, combined with the dip structure of the in- 
terconfiguration spectral function (J 1 — / 2 or p rg (u;)), leads to a bendover in resistivity at low 
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temperatures in the one-channel parameter regime within our simplified model. We believe 
this feature will go away in more realistic models. For example, another contribution to 
the Tg atomic spectral function comes from the convolution between f° and f 1 J = 5/2 T$. 
This spectral function will not be depleted at the fermi level, but instead will build up its 
spectral weight due to the weak Kondo resonance structure just above u = 0. Thus the low 
temperature bendover in the resistivity will disappear. 

The channel- dependence of the spectral functions, combined with the angular depen- 
dence, leads to a distinct behavior in the angle-averaged conduction electron lifetime for 
different numbers of channels. The spectral depletion becomes pronounced in the one- 
channel regime due to the angular average over the T 7 "hot spots." The development of 
the dip structure just below the Fermi level is correlated to the resistivity bendover at low 
temperatures. This feature does not occur in the two- and three-channel regimes. A Kondo 
resonance related peak develops near the Fermi energy, whose position with respect to the 
Fermi level depends on the channel character of the model parameters. Details are displayed 
in Fig. [T^. Note that the anisotropic scattering rate includes explicitly the hybridization as 
opposed to the other spectral functions like the atomic electron spectral function and the 
dynamic susceptibility. 

The transport coefficients are calculated using the Kubo formula under the assumption 
of dominant scattering in the I = 3 channel. Our results are defined in terms of transport 
integrals I n given by equations. 

df(oo) 



I n (T) = J du u n t(u,T) 



du 



(130) 



t(cu,T) = J^-rCk,co,T). (131) 
Here f(uj) is Fermi function. The resistivity is calculated using the equation: 

' — Jo(T). (132) 



p{T) m 

As we stated before, the resistivity shows a bendover at low temperatures in the one- 
channel regime within our simple model (not shown). This derives from the depleted spectral 
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weight in the f 1 — f 2 sector and zeros of Q^ 2 \k). 

In the two-channel and three-channel model parameter regimes, the resistivity initially 
increases logarithmically through T and then saturates with a power law to a constant 
with further decreasing temperatures. The resistivity near zero temperature obeys a scaling 
behavior as shown in Fig. [T7| confirming our zero temperature analysis with the scaling di- 
mensions, A 2 = 1/2, A3 = 2/5. These results agree with the conformal field theory results 
gg|, p(T) = p(0) [ 1 - a[T/T ] A ™ ] for T < 0.06T , A n = 2/(n + 2) for the overcom- 
pensated multi-channel Kondo models (note that the power law exponent is independent of 
the impurity spin size). We note that the region where strict T A " behavior holds is below 
about 0.05 To. A fit to the resistivity of Ce x Lai_ :r Cu2.2Si2 (see Fig. [[]) is good until low 
temperatures where the data breaks below theory. This suggests a possible crossover to a 
new fixed point which could be set by intersite interactions (producing a spin molecular 
field) or a weak noncubic symmetry for the Ce 3+ ions. 

The thermopower is a sensitive measure of the asymmetry in the scattering rate and 
the density of states (DOS) near the Fermi level. Since we are assuming the symmetrical 
Lorentzian/ Gaussian DOS for the conduction band, the sign of the thermopower is deter- 
mined by the asymmetrical scattering rate. 

^ = ~h m (133) 



As shown in Fig. [18], the thermopower shows explicitly the distinct behavior for different 
numbers of channels varying from the two-channel regime, through the three-channel , fi- 
nally to the one-channel regime. The hybridization between f° and f 1 leads to the Kondo 
resonance peak just above the Fermi level and more spectral weight above the Fermi level. 
Since particle scattering dominates, holes are the main carriers in the one-channel cases, 
thus leading to a positive thermopower. On the other hand, negative thermopower arises 
in the f l — f 2 sector. In this case, hole scatterings dominate and particles are the main 
carriers. In the three-channel regime, energy structures are symmetric. However, the double 
degeneracy in the f 2 configuration leads to a weak hole scattering dominance resulting in a 
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negative thermopower at low temperatures. The overall magnitude of the thermopower is 
slightly reduced when anisotropy at the cubic sites is included. 

VI. DISCUSSION AND CONCLUSION 

We have introduced and studied a realistic model Hamiltonian for Ce 3+ impurities with 
three configurations (/°, f 1 , f 2 ), which are embedded in cubic normal metals. This simple 
model shows competition between the Fermi liquid fixed point of the one-channel S — 1/2 
Kondo model and the non-Fermi liquid fixed point of the two-channel S — 1/2 Kondo model. 

We studied a simplified Anderson model using the NCA. This simple model covers one-, 
two-, and three-channel Kondo physics depending on the model parameters. All the calcu- 
lated physical quantities show the signatures of the Kondo effect appropriate to the differ- 
ent numbers of relevant channels. The magnetic susceptibility agrees with the exact Bethe 
Ansatz results for the two- and three-channel model parameter regime and has the correct 
scaling dimension agreeing with the conformal field theory results. Entropy and specific 
heat calculations also agree with the Bethe ansatz results within 0(1/N 2 ) approximation. 
Though the resistivity in the one-channel regime bends over with decreasing temperatures, 
we do not believe that this feature will survive when all the energy level structures, especially 
f 1 J = 5/2r 8 , are included. In the two- and three-channel parameter regimes, the resistivity 
increases logarithmically and saturates with decreasing temperatures. The low temperature 
behavior leads to power laws in agreement with the conformal field theory results. 

Since the thermopower is very sensitive to the density of states structure and the scatter- 
ing mechanism near the Fermi level, its sign and its magnitude depend on which fixed point 
is stable at zero temperature. In the one-channel regime, electrons are strongly scattered 
due to the Kondo resonance above the Fermi level and the thermopower remains positive 
definite and large due to the resonant scattering. In the two-channel case, holes are scat- 
tered off the impurity sites stronger than electrons. Thus electrons are the main carriers 
leading to the negative thermopower. In the three-channel regime, though hole scattering 
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is reduced compared to the two-channel case, electrons are still the main carriers due to the 
degeneracy imbalance between the singlet f° and the doublet f 2 J = 4r 3 . The thermopower 
remains negative with reduced value. We calculated the dynamic magnetic susceptibility 
and characterized it with its peak position (T(T)) as a function of temperature. We can 
see the clear difference between one-channel and the overcompensated cases. While T(T) 
decreases and saturates to a constant value of order Tk with decreasing temperatures in 
the one-channel case, T(T) goes to zero almost linearly with decreasing temperatures in the 
two- and three-channel case. 

We now discuss the experimental relevance of our model study to the Cei_a;La x Cu2.2Si2 
alloy |I7 |. We already stressed the experimental evidences supporting the two-channel Kondo 
effect in this alloy system in the section I. The thermopower for CeCu2Si2 changes sign 
around 70 K and is negative and large below |18],[45|. As our numerical calculation shows, 
the thermopower is negative and large in the two-channel regime. This result compares 
well with the experimental findings for the stoichiometric system with x — 1. We believe 
the sign change comes from the Kondo resonance of f° and / 1 r 8 , which lies above the 
Fermi level. Further experiments are required for the alloy system with excess Cu. For 
comparison, we note that CeAl 2 or CeAl 3 have a positive thermopower large compared 
to transition metals at high temperatures and have a sign change at low temperature which is 
still bigger than the Kondo temperature. Our thermopower calculation and the thermopower 



dependence on the unit cell volume |46| suggests that the alloy system Cei_ x Laa;Cu2.2Si2 can 
go through the three-channel model parameter regime with external pressure. Renormalized 
atom calculations further suggest a destabilization of f 2 relative to f° with initial increasing 
pressure 
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In addition to the thermopower, the future neutron scattering experiments for 
Cei_ a; La 3 ;Cu2.2Si2 alloy can search for the dependence of T(T) (the peak position of the 
dynamic magnetic susceptibility) consistent with our proposal for the channel number. 
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APPENDIX A: NCA SELF ENERGY IN THE PRESENCE OF MULTIPLE, 

SAME IRREPS 

In the f 2 configuration, there are 9 irreducible representations of T 3 in a cubic symmetry. 
NCA self energy diagrams should be generalized to get the right Kondo energy scale. As 
an example, we may consider the mixing process between f 1 J = 5/2T 7 and / 2 r 3 's. The 
relevant NCA integral equations are 

2 



EH = - £ yd^T" J de f(e)N(e) D %J {e + u), (Al) 

ij 

IL» = ^ yl^r" J de f(-e)N(e) G(e + u). (A2) 



The conduction electrons of only Tg symmetry are involved in the mixing process. The 
/ 2 r 3 's Green's function now becomes a 9 x 9 matrix. This generalization can be seen most 
clearly by looking at the self energy diagrams of / 2 IVs. The incoming T 3 does not need to 
be the same as the outgoing T 3 . This leads to the matrix Green's function for / 2 r 3 's. 
According to the Schrieffer- Wolff transformation, the effective Hamiltonian is 

ffi = 4,(0); J = j2^^. (A3) 

n i=l e 2 e l 

Now we can show that the above NCA integral equations lead to the right Kondo energy 
scale to leading order. 

G(u) -> -. (A4) 

v ' uj - e l + i6 

With this replacement, the T 3 self energy matrix becomes 



37 



Ilij(w) = -jTiTj log 

7T v 



UJ — €\ 



(A5) 



UJ — D — t\ 

for a symmetric flat conduction band DOS with a half width D. Now the Kondo energy 
scale is determined by 

det [ (u - ei)5ij - Ti^uj) ] = 0. (A6) 
With the substitution u = e\ — T , we can find the following 



log 



D + T Q 



iV(0)J' 



(AT) 



N(0)J=l±-fi-. (A8) 

71 i=l C 2 _ e l 



Here I\ = ttA^(O) | V^ 2 1 2 . 

In conclusion, we have shown that the inclusion of all the r 3 's in the f 2 configuration 
leads to the enhanced two-channel exchange coupling. 

APPENDIX B: ZERO TEMPERATURE ANALYSIS: OVERCOMPENSATED 
MULTICHANNEL ANDERSON MODEL. 

We discuss the following NCA integral equations at zero temperature. 

1 N T r ~ 

G g (z) = -— ; E g (z) = — / deN(-e)f(e)G x (z + e), (Bl) 

1 N T r 

GM = — -; ^(z) = -5- / deN(e)f(e)G g (z + e). (B2) 

z - e x - h x (z) 7T J 

Here N(e) is the normalized DOS such that N(0) = 1. For definiteness, we assume that 
the excited state has one less electron than the ground state. Before diving into the zero 
temperature analysis, we point out the applicability of NCA approach to the Anderson 
impurity model. Whenever the ground and the excited states have degeneracies of N g and 
iV x , respectively, the above form of NCA integral equations are obtained. The above NCA 
integral equations also derive from the A^-channel, Sj = (N g — l)/2 models. This is an 
artifact of the NCA approach. The NCA can not distinguish between these two different 
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models. Note that not all the Anderson models with the degeneracies of N g (the ground 
state) and iV x (the excited state) map into the iV x -channel, Sj = (N g — l)/2 models. The 
Schrieffer- Wolff transformation is essential to see this connection clearly. As an example, 
the excited triplets in the f 2 J = 4 in our Ce 3+ model leads to the one-channel exchange 
interaction of the impurity f l Y 7 pseudo spin Sj = 1/2 coupled to the S c = 3/2 conduction 
electrons [34f instead of three-channel Kondo model. With this restriction in mind, we now 



study the zero temperature analysis of the above NCA integral equations [0,^4 
The dynamical quantities we are interested in are 

p{w) = I de[ a x (e)A g (e + u) + A x (e)a g (e + u)], (B3) 



XgH = J de [ a g (e)A g (e A g (e)a g (e + u)], (B4) 

a g (u) = -—A g (u); a x (u) = -—A x (u). (B5) 
z f z f 

We can show that 

a g (u)\G g (iu)\- 2 = ^ j deN(-e)f(-e)a x (uj + e), (B6) 
a x (u)\G x (uj)\- 2 = [ deN(e)f(-e)a g (u; + e). (B7) 

7T J 

At zero temperature, Fermi function is reduced to step function. Thus, the self energy 
equations are simplified: 

N T r° 

E g (z) = ^ deN(-e)G x (z + e) } (B8) 

7T J-oo 

E JC (z) = - 5 -/ deN(e)G g (z + e), (B9) 

TT J-oc 



n x t 



oo 



a g (u)\G g (Lu)\- 2 = J de N(-e)a x (uj + e), (BIO) 
a x H|G x Mr 2 = / de N{e)a g (u + e). (Bll) 

TT JO 

For a flat conduction band with cut-off, [— D,D], we can reduce the above equations to 
differential equations: 

g g {u) = -1/G g (u); g x (u) = -1/G x (u), (B12) 
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d 

t 9 ^ = 



I -^iE-L-; g g (—D) — D + e g , 



5x(^) : 
NX 1 



a g(^)I^HP 

a x (w)|ftc(w)| 5 



iv x r 

7T 
7T 



7T 9 g (w)' 

a x (cj), 
a g (u). 



g x (-D) = D + e x , 



(B13) 

(B14) 

(B15) 
(B16) 



Removing the u dependence, we find the relationship between g g and g x , in terms of the 
integration constant which connects the low energy and high energy states. 



9g 



D + e s 



cxp 



iv g r 



cxp 



D + e^ 



iVx/TVg 



(B17) 



Since the zero temperature analysis is meaningful only when D ^> e gjX , we will replace 
D + e gjX by D. That is, 



9g 



= exp 









iV x /iV g 






.A. 





A" 


N x /N g 


"7r(e g -e x )" 




exp 





And 



T 
T = D 

A = £. 

7T 



[ JV g a g (c<;)0 g (c<;) + A?" x a x (^)5( x (cj) ] = N g a g (uj) + N x a x (uj). 



(B18) 
(B19) 
(B20) 



(B21) 



Due to the sharpness of Fermi function at the Fermi level, the spectral functions have a sharp 
cut-off at the threshold energy E . Since A n (uo) vanishes below this cut-off energy, T, n (u) is 
purely real and g n (ou) does not vanish or cross the frequency axis below the threshold energy 
E . That is, g n {uj) is positive definite below the threshold energy E . 



1. Leading asymptotic behavior 

We can derive the asymptotic behavior near the cut-off energy. Since a s (u) and g x {uj) 
vanish at u = E , we can approximate the above equations: 
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Eq — CO 
Eq — UJ 

T 



= / dy 
Jo 



9 s (y) 



9g (y) + N g A N g A Jo 
o V g x {y) + AT X A ~ A^ X A Jo 



dy g g (y), 
dy g x (y), 



A 



(B22) 
(B23) 
(B24) 



Here A = T/ir. From the above, we can find 



To 

?x(w) 



i^r g , 



T = £ 



AT 



iV g + iV x ' 



cxp 



7r(e g - e x ) 



iv g r 

N s + N x ' 



Q = (iVg + iV x ) 



E — uj 
To ' 



(B25) 
(B26) 
(B27) 
(B28) 
(B29) 



Since the zero temperature analysis is based upon the assumption |e n | <C TJ, the realistic 
Kondo energy scale is given by the replacement of D + e n — > TJ. The asymptotic behavior 
right above the cut-off £o can be found from the the expressions below E by the analytic 
continuation: 



g s (u + i5) 
To 

g x (uj + iS) 
A 

A g (w) 



e -ia s n 
1 



7tT 
1 



sin(a g vr) - £ ), 



A x (o>) « — sin(a x 7r) ax 6>(w - £7 ). 

7rA 



(B30) 

(B31) 
(B32) 
(B33) 



Here the phase was determined such that the spectral function is positive definite above the 
threshold energy. Since [ N g a g g g + N x a x g x }u)=e = 1, we can deduce that 

1 



a g (u) 



n- a z9(E -uj), 



(N s + Ay T 



(B34) 
(B35) 
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From the above asymptotic expressions, we find 

11 f sin(a g 7r) 9(u) 

r (iV g + Ax) 2 [ sin(a x 7r) 6{-u) 

= i [N g l Nx r (B36) 

** {u) * ^ (^+T X ^ B(a - ax) sgnM '" r2Qs ' (B37) 

u, = (N s + N x )^. (B38) 



Here q) is the beta function. 



2. Next leading asymptotic behavior 



From the above analysis, we can see that it is more appropriate to use the dimensionless 
quantities. Here we collect all the relevant formulas from the above. 



d_ 

d 



_d_9g 
dttT 

d g x 
dVL A 

a x (u)\g x (u))\ 

E — uj 

Eq—UJ 

To 



1 A 
Ax + ~g~ x 



N S A g s 
a g A T a x (u), 

a x A Tq a g (oo), 



9 *dy- 9 * {V) 



■ 9g (-D) = D + e g , 
; g x (-D) =D + e x , 



= exp 



o y g s (y) + N g A' 

o y g x (y) + N x A , 
ir(g g - g x 



Expanding the last relation, we get 



9g 
To 





7V x /AT g / 


.A. 





1 + 



Na-T 



To 





9* 




.A. 



#x 



AgA AgA 



A 



N x /N g 



+ 



Then it is straightforward to show that 

^ = f2 a « \l - 9l Q Qg + 2(? 2 + 
Jo L 



(B39) 
(B40) 
(B41) 
(B42) 
(B43) 
(B44) 
(B45) 



(B46) 



(B47) 
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1T X 



1 + 2^i !T g - # 2 1T X + 
1 



9i = 



T 



N g + 2N X A' ^ 2N g + N x 
Furthermore, writing in the Taylor expansion form, 

1 



a 



;i " [ N g + N x ] T 
1 



[ iV g + iV x ] A 
we can find the following relations 



ft -0 * 
fT Qx 



1 + ai fi Qg + a 2 1T X + • • 
1 + 6i fi ag + b 2 fi a * + • • 



9(Q), 

9(d). 



6i = 2ai - 45-1 ; a g 6 2 = a 2 + 4y 2 , 
a x ai = 6i + 45(1; a 2 = 26 2 - 4# 2 . 



Finally we find 



1 



N s + N x ] T 
1 



a x = 



fT Qx 



[N s + N x ] A 
And the pseudo spectral functions are 
1 



i + ofi Qg -4# 2 rr x + ---J 9(Q), 
1-4^1 Q Qg +0 fi Qx + 



' " n-V;, 

sin(a g 7r) 



sin(a g 7r) — 2g 2 sin[(a g — a x )ir] \Q\ a * + ■ ■ ■ 9(—il) 



7tT 

A ^ sin(a g 7r) 



VTA 



lfil- Qg 


1 + 4g 2 cos(a g 7r) 


|fi| ax 


+ ••• 




-Q), 


\tt\- a * 


1 — 4gi cos(o; g 7r) 


|fi| Qg 


+ ••• 




-f2). 



(B48) 
(B49) 



(B50) 
(B51) 



(B52) 
(B53) 



(B54) 
(B55) 



(B56) 
(B57) 



The scaling dimensions found here all agree with those found in the conformal field theory 
for the overcompensated cases. 



3. Physical quantities 

Using the above results, we can find the dynamic susceptibilities for the ground and the 
excited configurations 
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XM = sgn(cu) 



sm [a„n 



x 



(7V g + N X )HT 

B{a x , a x ) la)! 1-208 -4g 2 [l- cos(a g 7r) ] 5(a x , 2a x ) |cD| 2_3ois + 
T sin(a x 7r) 



(iV g + iV x ) 2 7rA 2 



x 



B(a s ,a s ) \Q 



~,|l-2a x 



4gi [ 1 — cos(a x 7r) ] B(a g ,2a g ) \lo 



~,|2-3a x 



+ 



(B58) 



(B59) 



These correlation functions are directly proportional to the imaginary part of the corre- 
sponding dynamic correlation functions. The functional forms are marginal Fermi liquid 
type for the case N g = N x , which correspond to the overcompensated case. The convoluted 
atomic electron spectral functions are 
sin(a g 7r) 



p t (w > 0) 



pi(u < 0) 



, (B60) 



(iV g + iV x ) 2 7rA 

x [B(a x ,a g ) — Agi B(a x ,2a g ) \Cu\ ae + Ag 2 cos(o;g7r) B(a g ,2a x 

sin(a; x 7r) 
(A^ g + iV x ) 2 7rA 

x [B(a x , a g ) + Agi cos(« x vr) B(a x , 2a g )\uj\ as - Ag 2 B(a g , 2a x ) \u\ a * H ] . (B61) 



One important observation is that the scaling dimensions agree with those obtained from the 
conformal field theory treatments for the overcompensated multichannel S = 1/2 models. 
From this result, we can argue that the convoluted local electron spectral function is peaked 
right at the Fermi level (Kondo resonance) for the overcompensated models. This result 
seems to be independent of the occupancy of the ground level. In fact, the numerical results 
with NCA confirm this conclusion. Hence we conclude that Kondo resonance peak in the 
overcompensated models sits right at the Fermi level irrespective of any model parameters. 
In the single channel Anderson model, the position of the Kondo resonance peak is adjusted 
by the Friedel sum rule (Fermi liquid ground state). That is, the occupancy of the ground 
configuration determines the Kondo resonance peak position. 

We can also generate the low temperature dependence of some physical quantities. Re- 
sistivity and thermopower can be evaluated using the Kubo formula. 



p(T) = p(0) 1 - c 



T 



; v = min(a g , a x ), 



(B62) 
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(B63) 

APPENDIX C: CONDUCTION ELECTRON SCATTERING RATE. 

Here we derive the conduction electron scattering time in the presence of the Anderson 
magnetic impurities. In general, the conduction electron scattering rate (the inverse of the 
conduction electron lifetime) is different from the transport scattering rate. This difference 
derives from the vertex correction when we calculate the current-current response function. 
The simplifying feature of the Anderson model is that only one partial wave state is coupled 
to the atomic electron state for an isotropic hybridization. In this simplified model, the 
transport scattering rate is proportional to the imaginary part of the conduction electron 
self energy. 

The conduction electron scattering time is given by the angular average of the anisotropic 
one. 

/dk a ^ a 
^— k^k u T^ u (k,uj). (CI) 

Since t m „(— k, u) = r^ v (k,uj) in the Anderson model, the scattering time matrix becomes 
diagonal. Thus we have 

/dk a ^ a 
— kykp r(k,w). (C2) 

In the dilute Anderson impurity limit, the conduction electron self energy can be approxi- 
mated by 

E(te, iuj) = N imp J2 I < ka\V\i > | 2 G^iuj), (C3) 

i 

T-^ka,^) = 27cN imp J2 \ < ka\V\i > | 2 pt(u). (C4) 

i 

Here the index i is the atomic electron's good quantum number. We can write the self 
energy as 



Q(T) oc 
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1. With LS coupling, but without CEF. 

Y,(ka, iuj) = N imp I < ka\V\jm > | 2 G t (J; iu), (C5) 

jm 

Gf(J; r) = -< T T f jm (r)fj m (0) >, (C6) 
<A;a|l / |jm> = ^ < ka\V\m 3 (3 > < 3m 3 ; ^/3|jm >, (C7) 

rrt 3 /3 

The total angular momentum is a good quantum number in this case. 

2. With both LS coupling and CEF. 

H(ha,iw) = N imp I < ka\V\r®d c > | 2 G t {Y^;iw), (C8) 

Gf(r?>; r) = - < T r / r , j)dc (r)/^ )dc (0) >, (C9) 

< ka\V\T®d c > = < ka\V\jm > < ]m\Yfd c > . (CIO) 

m 

The CEF irreps are good quantum numbers. 
In the above, two successive unitary transformations have been used. 

fm :i a = < 3m s; 7>a\jm > f jm (Cll) 

jm 

fjm = E < jrn\F®d c > f r u W (C12) 

cd c 

For the isotropic, spin-independent hybridization and the free electron conduction band, 
the mixing matrix is given by 

< ka\V\m 3 [3 > = ^ V(k)Y 3m3 (k) 5 a/3 (C13) 
V(k) = v 7 !^ (-if J dr r 2 j 3 (kr)V(r)R{(r). (C14) 

< ka\V\jm > = W— V(k) < ka\jm > (C15) 

< ka\jm > = Y, Y 3m 3 (k) < 3m 3 ; ]-a\jm > (C16) 

m 3 1 

< ka\V\Ti j) d c > = ^ V(k) < ka\r®d c > (C17) 

< ka\Tfd c > = £ Y 3m3 (k) < 3m 3 ; \0i\3m >< jm\V^d c > . (C18) 



m 3 m 
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Here f2 is the volume of the system. Hence the conduction electron self energy in the isotropic 
hybridization can be written as 

T,(ka,iu) =47m imp |V(A;)| 2 ^| < ka\jm > \ 2 G f (j;iu), (C19) 

jm 

= 4vm imp \V(k)\ 2 Y / \<ka\T^d c >\ 2 G { (T^;icu). (C20) 

jcd c 

Here only the diagonal elements of Green's function are nonvanishing. 

1. Without CEF 

When the CEF is neglected, the conduction electron self energy in a dilute impurity limit 

reads 

E(ka,iu) =47rn imp \V(k)\ 2 J2® U \k) G f (j;iu;), (C21) 

3 

T-\ka,u) = M 2 n- imp \V{k)\ 2 Y.® {l \k) PfO», (C22) 

j 

Q^\k) = ]T I < ka\jm > \ 2 . (C23) 

m 

Pf(j;uj) is the measurable spectral function for the atomic electrons with the total angular 
momentum j. 

The relevant angular functions defined above are, in fact, constants. 

e^(jfe) = ]T | < ka\jm > | 2 = ]T I < ka\Yfd^ > \ 2 = (C24) 

2. With CEF 

In the presence of the CEF, the CEF irreducible representations are good quantum num- 
bers. The atomic / electron operator can be decomposed into CEF irreducible representation 
components. Hence the conduction electron self energy in the dilute impurity limit is 

E{ka,iu) = 47m imp \V(k)\ 2 ]T Q^(k) G t (F^;iu), (C25) 
G$(k) = E I < ka\Yfd c > | 2 . (C26) 
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Here 0^ (k) is crystal harmonics. The anisotropic relaxation rate reads 

r-\ka,cu) = 8n 2 n imp \V(k)\ 2 ]T eg(fc) pt(lf;u). (C27) 

Here pi(T^;uj) is the measurable spectral function for the atomic / electron of sym- 
metry. This spectral function is given by the convolution of two neighboring configuration 
Green's functions. 

In the cubic crystal symmetry, the relevant crystal angular functions are 

e? /2) W = EMT/lir? /2) «>l 2 

a 

' [ 35cos 4 #-30cos 2 #-5 + 5sin 4 #cos4v? ], (C28) 



And 



32tt 

e? /2) W = El<^T/lir? /2) ;^>| 2 

net 

= -r- [35cos 4 #-30cos 2 # + 19 + 5sin 4 #cos4v?]. (C29) 

327T 



©6 Cr/2) (*) = EI<*T/i|r 6 cr/2) ;a>| 2 



7 34 
[ 50 cos 6 9 - 30 cos 4 9 - 10 cos 2 9 + "' 



256tt L 3 

-10 sin 4 #(5 cos 2 l)cos4^ ], (C30) 

ep/ 2 )(fc) = EK*T/i|r? /2) ;a>| 2 

'° sin 2 0[ - 14cos 4 # + 28cos 2 # + 2 



256tt 

-2sin 2 #(7cos 2 fl + l)cos4y? ], (C31) 

e 8 cr/2) (*) = EI<*T/i|r 8 cr/2) ;na>| 2 

net 

1 55 
' ' 70cos 6 # + 105 cos 4 - 40 cos 2 # + — 



32tt L " 3 

+5 sin 4 0(14 cos 2 fl- l)cos4y? ]. (C32) 

Note that the above crystal angular functions are normalized and satisfy the correct sum 
rules. 
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J dk e$(k) = 1; J dk e { i\k) = 2, (C33) 

EfW~- (C34) 

Furthermore, the crystal harmonics for j = 5/2 multiplets can be rewritten as 

Qf /2 \k) = ±- [ 6 - *(*) ]; 0? /2) (fc) = ^ [ 6 + *(*) ], (C35) 

$(fc) = 15 cos 4 0-10 cos 2 + 1 + 5 sin 4 cos 2 2<p. (C36) 
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FIGURES 



FIG. 1. Comparison of experimental resistivity with our numerical calculation. Our 

numerical calculation shows a T 1 / 2 behavior at low temperatures. Experimental resistivity (points 
with To = 10 K, from Ref. [lTj) is compared with our numerical results. Three different symbols 
mean a set of three model parameters (□ for model set 1; O f° r model set 2; A for model set 3). 
The low temperature deviation suggests a possible crossover to a new fixed point. 
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FIG. 2. Crystal electric field energy level scheme for f ,/ 1 ,/ 2 configurations. The 

one-channel & two-channel Anderson model Hamiltonians are developed from these CEF energy 
states. The one-channel Kondo model derives from the f° singlet and the magnetic f 1 T^ doublet 
which mix through hybridization with the T7 conduction electrons. The two-channel Kondo model 
derives from the magnetic ^T-j doublet and the nonmagnetic f 2 Ts doublet which mix through 
hybridization with the Tg conduction electrons. 
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FIG. 3. Schematic digrams of two channel degrees of freedom. The upper figure is for 
Tg+ and the lower figure for Tg— of J = 5/2 conduction electron partial waves. 
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FIG. 4. Scaling diagrams of one-channel & two-channel Anderson model up to third 
order. Dashed lines are for the Ce 3+ impurity Tf pseudo spins and the solid lines are for the 
or r§ n conduction electrons. The diagrams labeled as above lead to the scaling equation for the 
one-channel exchange coupling, J\. To obtain the scaling equations for the two-channel exchange 
coupling, J2, the labels T7 for the external conduction electrons (solid lines) should be replaced by 

T8n- 
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FIG. 5. Leading skeleton self energy diagrams from the NCA. Since two different 
symmetry conduction electrons are involved in the hybridizations of f° — f 1 and f 1 — f 2 , our NCA 
self energy diagrams become simplified. The diagram (a) is the self energy for the /° atomic state 
(wiggly line). The diagrams (b) are for the / 1 T7 atomic state (dashed line). The diagram (c) is 
for the / 2 r3 atomic state (dotted line). The solid line is the conduction electron propagator: the 
first two are for the Tj and the second two are for the Tg conduction electrons. 
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FIG. 6. Channel dependence of the Entropy. The NCA calculation of entropy clearly 
shows the right magnitude of residual entropy depending on the relevant channel numbers. Solid 
lines are Bethe-Ansatz curves. Note that the temperature scale is linear for M = 1. Referring to 
Table I: (a) M = 1 case - □ for model set 8; O for model set 7; A for model set 6. (b) M = 2 
case - □ for model set 1; O f° r model set 2; A for model set 3. (c) M = 3 case - □ for model set 
4; O f° r model set 5. 
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FIG. 7. Channel dependence of the Specific heat. Solid lines are Bethe-Ansatz results. 
The comparison to the Bethe Ansatz is complicated by the background derived from the in- 
ter-configuration peak. In the M = 1 cases, the Kondo temperatures are adjusted such that 
the NCA numerical results fall on the exact Bethe Ansatz one. Thoughout all the figures pre- 
sented in this paper the estimated Kondo temperatures in the one-channel models are used except 
for the magnetic susceptibility curves (see Fig. ||). Symbols have the same meaning as in Fig. [6|. 
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FIG. 8. The scaling behavior of the static magnetic susceptibility. The static magnetic 
susceptibility obeys a scaling behavior for each parameter regime leading to the M = 1,2,3 fixed 
points. The agreement with the Bethe ansatz results (solid lines) is good for the overscreened cases 
[M = 2,3). Tq = Tk/0.3 for M = 2, Tk from Ref. jl3|). For convenience, we multiply x(T) by 
2.0 for M = 2. M = 1 case - O for model set 8; * for model set 7; * for model set 6; M = 2 
case - □ for model set 1; O f° r m odel set 2; A for model set 3; M = 3 case — h for model set 4; 
x for model set 5. The ground state phase diagram for the model in exchange coupling constant 
parameter space is drawn in the inset, where g% = N(0)J{, N(0) being the conduction band density 
of states at the Fermi energy. The solid diagonal line is for M = 3. 



61 



FIG. 9. Atomic spectral functions in the one-channel regime, poi is the interconfigura- 
tion spectral function which is obtained from the convolution between f° and f^Ti states, pu is 
the interconfiguration spectral function which is obtained from the convolution between f 1 Tj and 
/ 2 r3 states. One-channel Kondo effect leads to the Kondo resonance development in poi just above 
the Fermi level and the spectral depletion in pu right at u = 0. Spectral functions are displayed 
for model set 8. The temperature variations are T/D = 3.678 x 1(T 2 , 1.077 x 1CT 2 , 3.155 x 1CT 3 , 
9.239 x 10~ 4 , 2.706 x 10" 4 , 7.924 x 10" 5 , 2.321 x 10" 5 . 



62 



FIG. 10. Atomic spectral functions in the two-channel regime. Two-channel Kondo 
effect leads to the Kondo resonance development in p±2 at the Fermi level (T = 0) and the spectral 
depletion in poi right at u = 0. Spectral functions are displayed for model set 1. The temperature 
variations are the same as in Fig. ||. 
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FIG. 11. Atomic spectral functions in the three-channel regime. For this parameter 
regime (model set 4), two spectral functions are equivalent in the asymptotic limit after a parti- 
cle-hole transformation. The temperature variations are the same as in Fig. [9[ 
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FIG. 12. Temperature dependence of the Kondo resonance peak. The Kondo resonance 
peak position {lok) shows a different temperature dependence on the relevant channel numbers, (a) 
One-channel case {pm): ojk decreases and saturates to a constant with decreasing temperature, 
(b) Two-channel case (P12) an d (c) Three-channel case (pi2) : decreases and tends to zero 
with decreasing temperature. Symbols have the same meaning as in Fig. |(| Note the different 
temperature ranges between M = 1 case (top) and M = 2, 3 cases (bottom). 
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FIG. 13. Total atomic spectral function. The high energy structure essentially does not 
depend on either the temperature or the relevant channel numbers. On the other hand, the 
relevant channel number dependence shows up in the magnitude of the Kondo resonance structure, 
(a) One-channel case (model set 8): the Kondo resonance peaks are cut off to be compared with 
the two- and three-channel cases, (b) Two-channel case (model set 1). (c) Three-channel case 
(model set 4). The temperature variations are the same as in Fig. ^. 
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FIG. 14. Variation of x (<^) with temperature. The reduced dynamic magnetic suscepti- 
bility x"(o>, T)/x"(T(T), T) is displayed and show a rough scaling behavior between two exrema. 
One-channel case: model set 8; two-channel case: model set 1; three-channel case: model set 4. 
The temperature variations are the same as in Fig. [9|. 
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FIG. 15. Peak position (T(T)) of the dynamic magnetic susceptibility. The temper- 
ature variation of T(T) depends on the relevant channel numbers. In the one-channel case (a), 
r(T) approaches a constant value as T — ► 0. In the two- or three-channel cases (b) and (c), 
r(T — ► 0) — > which is none other than Marginal Fermi liquid behavior. Symbols have the same 
meaning as in Fig. H Note the different temperature ranges between M = 1 case (top) and 
M = 2, 3 cases (bottom). 
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FIG. 16. Angle-averaged scattering rate, (a) One-channel case (model set 8): electron 
scattering is dominant. The cusp feature for M = 1 case right at the Fermi level derives from the 
spectral depletion of pi2 and zeros of the crystal harmonic @^ 2 \k) which is defined in the text, 
(b) Two-channel case (model set 1): hole scattering is dominant, (c) Three-channel case (model 
set 4): weak hole scattering dominance over the electron scattering comes from the degeneracy 
imbalance between the /° singlet and / 2 r3 doublet. The temperature variations are the same as 
in Fig. g. 
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FIG. 17. Low temperature dependence of resistivity. Correct power laws are found 
asymptotically for the finite temperature NCA calculation in the two- and three-channel cases. A 
power law of T 1 / 2 (T 2//5 ) is expected for M = 2(M = 3). Symbols have the same meaning as in 
Fig. | 
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FIG. 18. Channel dependence of the thermopower. The thermopower is a measure of 
asymmetry in the density of states and the scattering rate with respect to the Fermi level. Large 
thermopower derives from the Kondo resonance scattering. For M = 1, the holes are main carriers 
leading to positive thermopower. For M = 2,3, electrons are dominant carriers leading to a 
negative thermopower. The thermopower for M = 2 compares favorably with the experiments for 
the stoichimetric system, CeCu2Si2. For more details, see the text. 
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TABLES 



Model parameters for the Ce impurity. This set of model parameters covers the single, 
three, and two channel Kondo regimes. T = T i = r i2 is the hybridization strength for both 
f° — f 1 and f 1 — f 2 mixing, respectively. We scan channel number M = 1, 2, 3 according to 
whether e 2 > 0, < 0, = 0. The Kondo scale T is estimated as described in the text. 



Set 


M 


T/D 


ei/D 


e 2 /D 


k B T /D 


1 


2 


0.2 


-0.4 


-0.1 


1.9081 x 10 -3 


2 


2 


0.2 


-0.37 


-0.07 


1.5964 x 10 -3 


3 


2 


0.2 


-0.35 


-0.05 


1.3492 x 10~ 3 


4 


3 


0.2 


-0.3 


0.0 


1.5224 x 10~ 3 


5 


3 


0.2 


-0.4 


0.0 


6.9413 x 10~ 3 


6 


1 


0.2 


-0.3 


0.05 


1.1957 x 10~ 3 


7 


1 


0.2 


-0.3 


0.07 


1.6740 x 10~ 3 


8 


1 


0.2 


-0.3 


0.10 


2.3914 x 10~ 3 
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